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^ . The aim of the paper is to provide a constructive method for recov- 

' ering a unitary matrix from experimental data. Since there is a natural 

, immersion of unitary matrices within the set of double stochastic ones, the 

■ problem to solve is to find necessary and sufficient criteria that separate 

CO ' the two sets. A complete solution is provided for the 3-dimensional case, 

, accompanied by a test necessary for the reconstruction of a unitary 

• matrix from error affected data. 
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! 1 Introduction 
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An n X n matrix M is said to be double stochastic if its elements satisfy the 
relations 

n n 

'mij>0, ^ TOy = 1, ^m,y = l (1) 
i=i j=i 

All such matrices form a convex set called the Birkhoff's polytope The 
unistochastic matrices are a subset of the double stochastic ones defined by 

my = |[/„f (2) 

where [/ is a unitary matrix, i.e. satisfies the relation 

UU* ^U*U ^In (3) 

where * denotes the adjoint, and is the n-dimensional unit matrix. 
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It is well known that for n > 3 there are double stochastic matrices that are 
not unistochastic [2], and from a mathematical point of view there are a few 
interesting problems that deserve to be solved: 

i) Given a double stochastic matrix, M ~ {fTiij), i,j ~ l,' • • ,1^, what are 
the necessary and sufhcient conditions for M to be unistochastic. 

ii) Supposing M unistochastic, to what extent the matrix U is determined 
by M, i.e. how many solutions one could get. 

Ill) If M is unistochastic, how one can reconstruct the unitary U from the 
given data. 

The problem i) was completely solved only for the 3-dimensional case 
A characterization of the subset of double stochastic matrices that come from 
unistochastic ones, Eq.Q, for n = 3, was given in [3], 0]- For n > 4 only partial 
results are known, see e.g. As a matter of fact the theoretical physicists 

working on this problem were not aware of the embedding of the unitary ma- 
trices into the double stochastic ones, so in general they were not interested in 
obtaining necessary and sufficient conditions for the separation of the two sets. 
In this respect see [S], Introduction, where it is said: We are not concerned here 
with the consistency problem, which amounts to obtaining necessary and suffi- 
cient conditions on the set of numbers \Ujk\ for this set to represent the moduli 
of a unitary matrix. In this paper we will provide such necessary and sufficient 
conditions on the independent parameters entering the unitary matrix, when 
they are expressed in terms of the entries of a double stochastic matrix. We 
mention also that the Theorem 1 in the paper which provides necessary 
and sufficient conditions in the ti = 3 case, is only an existence theorem, and it 
does not provide a constructive method for recovering a unitary matrix from a 
double stochastic one. 

Concerning problem ii) it was shown that the generic situation for n > 4 is 
the existence of a continuum of solutions, see 0-111 and ^J|-C21- In particular 
complex Hadamard matrices, i.e. unitary matrices with equal moduli, \Uij\ = 
l/-y/ri, have been found that depend on arbitrary phases, see |13|-|17|. 

The solution of the mathematical problem m) for n = 3 can be easily ob- 
tained, as we show in the following, but its applications in high energy physics 
put some challenges since the measured numbers are affected by errors. Usually 
one starts with a theoretical model formalized by a unitary matrix, and the 
problem is to recover the unitary matrix from measurements of its moduli, or 
of some "angular looking objects", Uij Uik Uik Uij, where bar means complex 
conjugation. When |C/yP are measured one gets a numerical matrix as 

V = {V^j) (4) 

where every Vij is affected by errors, so we do not know how far from \Uij\ are 
the measured values Vij , and by consequence the double stochasticity relations 
([T)l are only approximately satisfied, if any. More generally, that means that 
we have to test the full compatibility between the data and the unistochastic 
property. Second, by supposing the data are compatible to the existence of a 
unitary matrix we need a reliable algorithm for an explicit recovery of a unitary 
matrix from data, because some parameters entering the unitary matrix may 
have a physical significance, and we want to know them, an example being the 
CP-violating phase from the Cabibbo-Kobayashi-Maskawa (CKM) matrix. 

Indeed the recovery of a unitary matrix from experimental data is a central 
problem in the electro- weak interactions jJSl where the (assumed) unitary CKM 
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matrix |19j plays a fundamental role. Hence the recovery problem of a unitary 
matrix from experimental data is not an academic problem, it has many prac- 
tical consequences. Let remind that there are two big Collaborations, BaBar in 
US, and Belle in Japan whose efforts are to measure as exact as possible the 
b-quark related entries of the CKM matrix. In Europe, at CERN, there is under 
construction the LHC machine, one of its main aims being a better understand- 
ing of the so called B-physics. Thus a reliable algorithm for recovery of unitary 
matrices, able to obtain the independent parameters of the CKM matrix, could 
have consequences on both the design of future experiments, as well as on the 
design of future high energy machines, including neutrino factories. 

The main goal of the paper is to provide a reliable algorithm for the re- 
construction of a unitary matrix from experimental data, when these ones are 
compatible with the theoretical model. We mention that the nowadays algo- 
rithm for reconstruction of a 3 x 3 unitary matrix from experimental data does 
not make use of the double stochasticity relations Q), and the phenomenolog- 
ical model used to describe the data and to reconstruct a unitary matrix from 
them is mainly based on the use of a single orthogonality relation, expressed 
as a triangle in the complex plane 1201-1101 j although the proof of Theorem 1 
from paper [H explicitly stresses the necessity of using at least two orthogonality 
relations. 

The paper is organized as follows. In sect. 2, starting from the properties of 
double stochastic matrices and the embedding relation Q of unitary matrices 
within that set, we describe the gauge group of unitary matrices, i.e. the group 
of the most elementary transformations whose action on the unitary matrices 
does not change the unitarity property, or the physical content. In sect. 3 we 
provide a parameterization of unitary matrices that is essential in devising a 
reconstruction algorithm in terms of physically relevant quantities. In sect. 4 we 
define two phenomenological models, unitarity condition method, and unitar- 
ity triangles method, by using the embedding ((SJ and the double stochasticity 
properties Wc find the necessary and sufficient conditions the data have to 
satisfy in each model in order to be consistent with the unitarity properties, and 
give the reconstruction algorithm of unitary matrices from double stochastic ma- 
trices. In sect. 5 we show that both the approaches are completely equivalent, 
then and only then, when they are formulated in terms of four independent 
moduli; a consequence will be that the second model has to use at least two 
orthogonality relations. In sect. 6 we describe the reconstruction algorithm of 
unitary matrices from experimental data that are compatible with the double 
stochasticity property. With this aim in view we define x^-tests that allow the 
recovery of unitary matrices from error affected data, including a method for 
doing statistics on moduli of unitary matrices. The paper ends by Conclusion. 



2 Unitary matrices and their gauge subgroup 



It is well known that an n x n unitary matrix depends on parameters, |31j , 
|32j ■ that are usually taken as n{n— l)/2 angles and n{n+ l)/2 phases, each set 
taking values within [0, 7r/2], and, respectively, [0, 2 tt). Eq.(|21) tell us that, given 
a definite unitary matrix, all the unitary matrices obtained by multiplying it at 
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left and/or at right by diagonal phase matrices, D = diagonal (6**^1, •• • , e**^"), 
with arbitrary real (f>i, i = 1 , • • • , n, generate a single double stochastic matrix. 
That means that we can simplify a little bit the form of a unitary matrix since 
the values of 2n — 1 phases are at our disposal, and a common choice for them is 
and/or tt. In high energy physics this property is known as phase invariance 
Thus we can take the entries from the first row and the first column 
as nonnegative quantities, such that the number of independent parameters en- 
tering a unitary matrix gets — (2n — 1) = (n — 1)^, and it is equal with the 
number of independent parameters entering a double stochastic matrix. Hence 
we could say that the embedding (O suggests that the "natural" coordinates to 
parametrize a unitary matrix could be the moduli of its entries. Unfortunately 
that can be done only for 3x3 matrices, but even in this case there are supple- 
mentary relations that have to be fulfilled by the moduli in order to get from 
them a unitary matrix. 

Besides these transformations there are another transformations: multipli- 
cation at left and/or right by permutation matrices. Permutation matrices are 
matrices whose elements on each row and each column are zero, but one that 
equals unity. They interchange rows, and, respectively, columns between them- 
selves. Both diagonal phase matrices and permutation matrices are subgroups 
of unitary matrices. If D denotes a diagonal phase matrix and P a permutation 
matrix then 

D D* = P P* = In 

Other equivalent unitary matrices can be obtained by taking the complex conju- 
gate matrix, and/or the transpose of the original one. If we denote the transpose 
operator by T, and by C the complex conjugation, these both transforms form 
a subgroup because 

= = Identity 

Thus the product group 

K^DxPxTxC 

is the gauge invariance subgroup of unitary matrices, and, in the following, we 
work only with the cosct defined by 

X = U{n)/K (5) 

where U{n) denotes the n-dimensional unitary group. The C invariance has 
an important consequence: the range of all independent phases entering U (n) 
is [0,7r]. The multiplicity of solutions appearing in the recovery problem of 
a unitary matrix from a double stochastic one will be given modulo the above 
simplest transformations of unitary matrices. 

3 Parametrization of unitary matrices 

For devising a recovery algorithm of unitary matrices from double stochastic 
ones, or from experimental data, we need an explicit parametrization of them. 
That will lead us easily to the necessary and sufficient conditions a double 
stochastic matrix has to satisfy for being also unistochastic, even if a complete 
explicit solution is found only for n = 3. These conditions will lead to separation 
criteria between the double stochastic and unistochastic matrices. 
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There are essentially two types of parameterizations: first the classical result 
by Murnagham, that states that the matrices from the unitary group U{n) 
are products of a diagonal phase matrix containing n phases, and n{n — l)/2 
matrices whose main building block has the form 

/ cos6i -sin6ie-*'^ \ , . 

\ sinee'^fi cos9 J ^ ' 

i.e. a paramctrization in terms of n(n — l)/2 angles 9i, and n{n + l)/2 phases 
ifii. The usual paramctrization of the CKM matrix |37j used in high energy 
physics is of this type. A second paramctrization, also through factorization, is 
that given in paper |38| . The idea behind such a paramctrization comes from 
the following sequence 

U(n) - X X ... X ^ X Uil) 

^ S^"-! X S'2"-3 X . . . X 5^ X 

sequence that shows that each factor can be parameterized by an arbitrary 
point on the corresponding complex sphere, i.e. by a single complex (n — k)- 
dimcnsional unit vector. Such a paramctrization could be appealing for high 
energy physicists since it shows that the information brought by each generation 
of quarks, or leptons is contained in a single complex unit vector. The first 
paramctrization of CKM matrix, , is of this form. The explicit realization 
of this paramctrization, that is the main result in |38j . is given by 

Theorem 1. Any element Un G U{n) can be factored into an ordered product 
of n matrices of the following form 

Un = B" ■ ■ ■ BI^_^ . . . B"~^ (8) 

where 



_( h \ 
^-'^ \ Bn-k ) 

and Bn-k G U{n — k) are special unitary matrices, each one generated by a 
single complex {n — k)- dimensional unit vector, bn-k G S'^'""'^-'"^. For example 
Bi = e}^ , where ip is an arbitrary phase. 

If Um G 5''^'"~^, m = 1, . . . , n, is parameterized by 

= (e"^i cos 6*1, e''^2 sin Oi cos 6*2, ... , e"^" sin Oi . . . sin Or^-iY 

where t means transpose, then the m columns of Bm are given by 



vi = y,n 



I e*"^! cos 6*1 

e*'^= sin 6*1 cos 6*2 



9i . . .sin Om-i j 
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and 

^k+i = -JTT viidi = ■■■ = 9k-i = 7r/2), A: = 1, • • • ,to - 1, 
at/k 

where in the above formula one calculates first the derivative and afterwords the 
restriction to 7r/2. 

In what follows we need an explicit parametrization of Bn-k, k = 0, . . . , n— 1, 
and by taking into account that wc work with representatives from the X coset, 
Eq.Q, we choose the corresponding n generating vectors as follows 

Un = (cos oi, sin «! cos a2, ■ ■ ■ , sin ai ■ ■ ■ sin a„_i)* 

y,i-i = (— cos 6i, e*'^i sin bi cos 62, ■ ' ■ 1 e^^"--^ sin 61 • • • sin 6„_2)* 

; (9) 

y2 = (—cos zi,e^^^ sin zi)* 
= -1 

With this choice the first row of [/„ has the form 

Mil = cos ai, U12 — cos bi sin ai, • • • ui„ ~ sin ai . . . sin zi (10) 

and the first column is given by ?/„. In the following wc will assume that C/„, 
Eq.®, has its first column given by ?/„, and its first row given by H10|l . and this 
is our standard form in the rest of the paper for a unitary matrix. A similar 
parametrization was obtained recently, see |39| : for the case n = 3 see also the 
papers [101 and gT]. 

4 Phenomenological models 

By "phenomenological model" we will understand in the following a relationship 
between the entries of a double stochastic matrix and the entries of a unitary 
matrix, the main goal being the finding of the necessary and sufficient conditions 
that separate the two sets. Depending on the context the "experimental data" 
will denote either the entries of a double stochastic matrix, or the numbers, 
affected by errors, measured in an experiment. Usually the experimental data 
on the 3x3 CKM matrix entries from the quark sector are given in terms of 
moduli of the unitary matrix that define the theoretical model, \Uij\, or angular 
looking objects UajUpkU akU pj, that are equivalent with the angles of the 
triangles generated by the orthogonality relations, where bar means complex 
conjugation. For the beginning we assume that the data have no errors, which 
is the current mathematical setting, i.e. the moduli are the entries of a double 
stochastic matrix, and we want to solve the problems i) — ni) from Introduction. 
For doing that we have at our disposal an explicit parametrization of unitary 
matrices, ©-01 a-nd the unitarity property 

UU* ^U*U ^ In 
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This concise form is equivalent with 2?! relations 



i=l 
i=n 

Y,\U^,?~l = 0. J = l,---,ri (11) 

i=l 

showing that the numbers = |C/ijP define a double stochastic matrix, that 
implies that only 2n — 1 relations from the set (|ll|l are independent, and by 
n{n — 1) orthogonality relations 

i—n 

Y^Uj^U;, = 0, j<k, i = l,---,n 

4=1 

i—n 

^t/„[/*,. = 0, j<k, z=l,---,n (12) 

i=l 

that can be visualized as polygons in the complex plane. The last relations are 
the supplementary relations the numbers Uij have to satisfy in order that the 
corresponding matrix should be unitary. The number of relations 1111(1 and l(12() 
is greater than n^, but we have written all of them since they could be useful 
in over-constraining the experimental data, that usually are affected by errors. 

The relations (|11() and, respectively, H12|l can be used to define two different 
phenomenological models. The first model is given by the relations l)ll|l together 
with 

m^j = \Uij\'^, i,j = !,■■■ ,n-'l (13) 

where rriij are the entries of a double stochastic matrix, and Uij arc the entries of 
a unitary matrix parametrized as in Theorem 1. The last relation is equivalent 
with the following equations 

mil = cos'^ fli, TOi2 = sin^ ai cos^ 6i, • • • , rnin = sin^ ai ■ ■ ■ sin^ zi, m2i = cos^ 02 sin^ oi, 

2 ■ 2 ■ 2 -2 -2 ("^^) 

m^i = cos a3 sm ai sm 02, • • • , m„„_i = sm cii ■ ■ ■ sm a„_i. 



771,22 = cos^ fli cos^ a2 cos^ 61 + cos^ 62 sin^ a2 sin^ bi + 
2 cos ai cos 02 cos bi cos 62 sin 02 sin hi cos /?i , 



(15) 



"T-32 — cos^ fli cos^ 03 cos^ bi sin^ 02 + cos^ 02 cos^ 03 cos^ 62 sin^ 61 + 

sin^ 03 sin^ 61 sin^ b2 — 2 cos ai cos 02 cos^ 03 cos bi cos 62 sin 02 sin 61 cos /3i + 

2 cos fli cos 03 cos bi sin 02 sin 03 sin bi sin &2 cos /32~ 

2 cos 02 cos 03 cos &2 sin 03 sin^ 61 sin 62 cos(/3i — /32 ) 7 etc. 

(16) 

where we have written only the simplest equations. It is easily seen from the 
above equations that, since mu and mn, i = 1, • • • , n— 1, are entries of a double 
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stochastic matrix, there is a unique solution for cosa^ G (0, 1), i = 1, • • • , n — 1, 
of the form 



COS ai = TOii, cos Ofc = ^— , fc = 2, • • • , 71 — 1 (17) 

1 - E,=i "la 

and similarly for cos &i , cos Ci , • • • , cos zi . Hence the number of angles that have 
to be found is [n - l){n - 2)/2 - (2n - 3) = (n - 2)(n - 3)/2. In that way 
the number of equations of the form (|15() - (|16|l we have to solve is only (ti — 2)^. 
We substitute the forms for cosa^, sina^, and those similar to, from the first 
column and the first row, in Eas. ll5|l - H16() . such that we get (ti — 2)^ equations 
that depend on (n — 1)^ moduli niy, z,j = l,...,n— 1. We do now a relabeling 
of the angles, hi — > i = 2, • • • , n — 2, C2 — *■ 6n-i, ■ ■ • , etc., and similarly 

for the phases. With this notation the necessary and sufficient conditions for a 
double stochastic matrix to be also unistochastic arc given by the relations 

(?i-2)(n-3) , , 

0<cos6, <1, z = l,--- ,^ 1 (18) 

(n-l)(ri-2) , , 

-l<cosA<l, i=l,--- ,^ '-^ (19) 

where cos bi and cos f3i are the solutions in terms of rnij of the (ri — 2)^ equations 
of the form (|15|I - H16|1 . The above relations are in the same time the separation 
criteria between the double stochastic and unistochastic matrices, and their 
fulfillment is equivalent with the existence of at least one unitary matrix com- 
patible with the moduli my. To check them we have to solve analytically or 
numerically the (n — 2)^ equations (|15|I - H16|I . Numerically this can be done 
when TTiij are the elements of a double stochastic matrix, but for the real case 
of experimental data with errors we need an analytic solution to be used in a 
X^-test, and until now this was found only for the case n = 3. In conclusion 
for 3x3 data coming from an exact double stochastic matrix we have only one 
constraint, namely, — 1 < cos/3i < 1- If the data come from an experiment we 
have to check also the compatibility of the entries from the first row and first 
column with the relations p7|l . i.e. to see if the conditions < cosa^ < 1 are 
satisfied. 

Taking into account the relation H13|l . rriij = \Uij\'^ , that shows how the 
unitary matrices are embedded into the double stochastic ones, the parameter- 
ization of unitary matrices by their moduli seems to be very appealing in this 
case, although it is not a natural one in the general case. A natural parameteri- 
zation would be one whose parameters are free, i.e. there are no supplementary 
constraints upon them, as Eas. H18|) - (|19|l . to enforce unitarity. 

The problem we rose in |11| was to what extent the knowledge of the moduli, 
uiij = \Uij\, of an 71 X n unitary matrix U determines U. If we identify the 
parameters to the moduli, they will be lying within the simple domain 

l? = (0,l)x...x(0,l) = (0,l)("-^'' 

where the above notation means that the number of factors entering the topo- 
logical product is (n— 1)^. We excluded only the extremities of each interval, i.e. 
the points and 1 that is a zero measure set within U{n) and has no relevance 
to the problem of recovery of a unitary matrix from a double stochastic one. 
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Nothing remains but to check if the new parameterization is one-to-one. A 
solution to the last problem is the following: start with a one-to-one parame- 
terization of U{n), as that given in the preceding section, and then change the 
coordinates taking as new coordinates the moduli of the (n — 2)^ entries; these 
ones are obtained by deleting the first and the last row, respectively, the first 
and the last column. The moduli of the first row and the first column are in 
one-to-one correspondence with the parameters entering the unitary matrix, see 
e.g. Eas. (|f 7|l . and the moduli entering the last row and the last column are 
uniquely determined by the double stochasticity property. Afterwords use the 
implicit function theorem to find the points where the new parameterization 
fails to be one-to-one. The corresponding variety upon which the application is 
not a bijective one is given by setting to zero the Jacobian of the transformation, 
i.e. 

J ^ 9(77122, ■ ■ . ,m2n-l, . . . ,m„_i„_i) ^ ^ ^^^^ 

d(bi, . . . , 6(„_2)(n-3)/2, Ply ■ ■ , /3(n-l)(n-2)/2) 

One gets that, generically, for 77 > 4 the unitary group U{ti) cannot be fully 
parameterized by the moduli of its entries, [H]-|17j. i.e. for a given set of moduli 
there exists a continuum of solutions, the simplest example being the case of 
complex Hadamard matrices ^SJ-EI- The maximum dimension of the above 
variety is (n— 2)^ — 1 = (n— 3)(r7 — 1). For 77 = 3, J 7^ 0, and only in this case the 
parametrization of a unitary matrix through the moduli could be one-to-one. If 
the moduli are outside of the above variety an upper bound for the multiplicity 

n(n-3) 

is 2 2 ^ bound that is saturated for 77 = 3, when there is essentially only 
one complex matrix, if we take into account the gauge invariance of unitary 
matrices. 

In the case of exact double stochastic matrices, as we showed before, only 
(77 — 2)^ moduli enter the game since the angles entering the first column and the 
first row are uniquely determined. To have a flavor of the problem we consider 
more in detail the case rt = 4, when there are four equations, two of them being 
Eas. H15|l - H16|l . and the last two are 

777 2 3 — cos^ 61 cos^ &2 cos^ ci sin^ 02 -I- cos^ fli cos^ 02 cos^ ci sin^ bi + 
sin^ a2 sin^ 62 sin^ ci — 2 cos ai cos 02 cos 61 cos 62 cos ci sin 02 sin bi cos /3i + 
2 cos 61 cos 62 cos ci sin^ 02 sin &2 sin ci cos 71 — 

2 cosfli cos a2 cos ci sin 02 sin 61 sin 62 sinci cos(/3i -\- 71) (21) 
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™33 = cos^ a2 cos^ as cos^ 61 cos^ 62 cos^ ci + cos^ ai cos^ 03 cos^ ci sin^ 02 sin^ 61 
cos^ 61 cos^ ci sin^ 03 sin^ 62 + cos^ 62 sin^ 03 sin^ ci + cos^ 02 cos^ 03 sin^ &2 sin^ ci 
+2 cos oi cos 02 cos^ 03 cos &i cos &2 cos^ ci sin 02 sin 61 cos j3i 
— 2 cos 02 cos 03 cos^ 61 cos 62 cos^ ci sin 03 sin 62 cos (/3i — /32 ) 
—2 cos ai cos as cos 61 cos^ ci sin a2 sin 03 sin hi sin 62 cos /32 
+2 cos a2 cos 03 cos &i cos^ 62 cos ci sin as sin ci cos(/3i — /32 ~ 7i ) 
+2 cos ai cos as cos 62 cos ci sin 02 sin as sin 61 sin ci cos(/32 + 71 ) 
+2 cos^ a2 cos^ as cos 61 cos 62 cos ci sin 62 sin ci cos 71 
—2 cos 61 cos 62 cos ci sin^ as sin 62 sin ci cos 71 
+2 cos ai cos 02 cos^ as cos ci sin a2 sin &i sin 62 sin ci cos(/3i + 71 ) 
—2 cosa2 cos as cos hi cos ci sin as sin^ 62 sinci cos(/3i — /32 + 7i) 
+2 cos a2 cos as cos 62 sin as sin 62 sin^ ci cos(/3i — /32) 
From Eas. (|14|) we get 



cosoi ~ J mil, sinai = vT— mn 



77191 /I ^ "ill — fn2i 
cosa2 = (/- = , sina2 = 1/ 

i — 77711 V 1 ^ "^11 



/ "Isi . / 1 - 77711 - 777-21 - 77731 

cosas = ,/- , smas 1/ (23) 

" 1 - 77711 - m-21 V 1 - "ill - m2i 



, , mi2 . , /I - 77711 - mi2 

cosbi = J- , sin6i = i/ 

i — 77711 V 1 ~ "^11 



777-13 . /I - 77711 - 777-12 - 77713 

cosci = \ - , smci = W 

1 - 77711 - m-12 V 1 - "ill - ™12 

By substituting relations H23() into equations (|15(l - (|16|l and H21|l - H22() we ob- 
tain four equations that depend on nine moduli 777ij, 7,j = 1,2,3, and on the 
non-relabeled parameters 62, /3i, P2, 7i 

/22 - m22 = 0, /23 - 77723 = 0, 

/32 - ?"32 = 0, /ss - 77733 = (24) 

In general the rank of the Jacobian matrix 

J _ d{f22, /23, /32, /33) 

a(&2,/3i,/?2,7i) ^ ' 

will be less than (77 — 2)^ since we know that there are particular solutions that 
depend on an arbitrary phase. If we look at the Jacobian (|25|l as a function 
on the Birkoff's polytopc, i.e. depending on (77 — 1)^ = 9 independent moduli 
r77ij it could be possible to find a domain where rank(J) = 4, i.e. in this case 
we have only one solution. Then the compatibility relations between the double 
stochastic matrix entries {rriij)^ j^^^, and the unitarity property, or in other 
words, the separation criteria between the two sets, are four, and they have the 
form 

0<cos&2<l, -l<cos/3i<l, -1< COS/32 <1, -l<cos7i<l (26) 



10 



The above relations are the necessary and sufhcient conditions the moduH 
of a 4 X 4 double stochastic matrix have to satisfy in order to exist a unitary 
matrix whose moduli coincide with rriij , and their intersection gives the maximal 
domain within the rriij simplex that is compatible to the existence of unitary 
matrices. 

In case when for a given numerical matrix. rank( J) < 4, there is no one-to- 
one correspondence between the entries ruij of a double stochastic matrix and 
the independents parameters entering a unitary matrix, i.e. there is at least one 
solution that depends on an arbitrary parameter, phase or angle. 

The checking of criteria 1)26(1 requires explicit analytic solutions for 

C0S&2, C0S/3i, COS/32, COS71 

in terms of rriij ^ and when these ones are not numbers, solving the Eas. (|24|l . 
is not a simple problem even with the symbolic calculation software packages 
nowadays available. The only results in this direction are those obtained in [Sj, 
however they have to be used with caution since the authors assumed that no 
matter how the numbers entering a double stochastic matrix are, the Eas. ((24(l 
have a physical solution. 

By taking into account the above considerations the following result holds 

Theorem 2. Suppose we have a generalized spherical coordinate system on 
the unitary group U{n), and let U € U{n) be a given matrix parameterized as 
in Theorem 1, through n{n — l)/2 angles, each one taking values in [0, 7r/2], 
and (n — l){n — 2)/2 phases taking values in [0, tt], and let M = (jUij) be 
a n X n double stochastic matrix, whose entries are supposed to come from a 
unistochastic matrix U by the embedding 

■niij = |C/,jP, i,j = !,■■■ ,n~l (27) 

From equations {14\ ) we get a unique solution for the angles entering the first 
column and the first row of U as follows 

2 2 '^21 2 

cos ai = mil, cos 02 = , • • • , cos a„_i = ^„_i , 

i - 1^11 1 - mil - Li=2 

cos 61 = ,---,cos zi = (28) 

i^-mii 1 - "^11 - 2]i=2 

We substitute them in Eqs. \21^ obtaining a set of in — 2)^ equations that, after 
relabeling of the angles and phases, is of the form 

fijibl, . . . ,b („-2K„~3) ,/?l,- • • , 13 (r;^-l)(n^) ^ rUi-j , 1 , j = 2 , ' • ' , 71 - 1 (29) 

The solutions of the above equations are compatible with the existence of a uni- 
tary matrix, if and only if, all the angles and phases satisfy the unitarity con- 
straints \lf^) - J^19\) . For n > 4 the solutions of equations I2fj\} could depend on 
arbitrary angles and/or phases on the variety obtained by setting to zero the 
determinant of Jacobian matrix of the transformation l^2y\) 

J _ d{f22, ■ ■ ■ , /2^^-l, /^^-12 ' ' ' , fn-ln-l) ^g^^ 

9(6i , • • • , b („-2)(n-3) , /3i , • • • , /3 („^i)(„-2) ) 
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Ifp is the rank of the Jacobian matrix the solution of \2y\) depends upon 
(n — 2)^ — p arbitrary parameters, angles and/or phases. Outside this variety 
the number of discrete solutions Ng satisfies 1 < Ng < 2 2 

Proof. Since wc use a spherical coordinate system the equations (|29() are 
trigonometric equations in our parameterization, as the example of case 77. = 4 
shows, Eas. (jl5(I . H16(l . (|2H) . H22(l . and consequently the multiplicity of the solutions 
may arise from the two possible phase solutions for all values of sine or cosine 
functions that satisfy Eas. 129(1 and the constraints (|18|I -H19 |I . The number of 
independent phases is (n — l)(n — 2)/2 and, since we do not make any distinction 
between U and its complex conjugate [/, condition that halves the number of 
solutions, the above bound for Ng follows. 

For n = 3 the Jacobian does not vanish and one gets 1 < A^s < 1, and this 
bound implies the existence of a complex unitary matrix, then and only then, 
when the values m^j , coming from a double stochastic matrix satisfy the relation 
-1 < cos/3i < 1. □ 

An example of a unitary matrix that cannot be recovered from its moduli 
is the following. If P and Qi, i = 1, • • • , m, are m x m and, respectively, n x n 
unitary matrices whose first rows and first columns are positive numbers and 
depend on p and respectively qt arbitrary phases, then the following mxn array 



/ PiiQi ■ ■ 



M 



(31) 



\ Plm Ql ' ' Pmm Qm j 

defines a unitary matrix that could depend on 



p + 91 + (m - 1) ^ q, (32) 



arbitrary phases. 

Indeed it is easily seen that we can multiply at left all the matrices (52, • ■ ■ 7 Qm 
by diagonal phase matrices -Dj — Diagonal{l, e'"^^'^ , ■ ■ ■ , , e"^"-i'J ), j = 2, • • • ,m, 
by preserving the entries from the first row and column of M positive numbers, 
obtaining a set of unitary matrices that are all applied in the same double 
stochastic matrix. Thus their recovery from the moduli cannot be done because 
all the arbitrary phases tpij disappear when one computes their moduli, and by 
consequence they do not appear in equations as ((T5|l - lfTB|) . 

A second phenomenological model can be defined by starting from the or- 
thogonality relations H12|l . but since in this case the polygons angles enter the 
game its formulation for arbitrary n is more difficult. Thus in the following we 
will discuss the case n = 3 for both the models, that has applications in high 
energy physics. 



4.1 Unitarity condition method 

In the n = 3 case the unitary (CKM) matrix is parametrized by four inde- 
pendent parameters given by the so called mixing angles, ^12,^13,^23, and the 
CP-violating phase S. Hence in the following we will change the notation from 
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section 3 to another notation that is more famihar to experimenters and phe- 
nomenologists. That means that in equations lO defining the generating vectors 
we make the substitution: (oi — > 6*12, 02 O23, bi 9i3, Pi ^ S); after that 
we make the notation 

cos dij — Cij, sinOij ~ Sij, ij — 12, 13, 23 

and by using Theorem 1 we get the following form 

C13S12 S12S13 
C23S12 -C12C13C23 - e"^si3S23 -C12C23S13 + e"'ci3S23 I (33) 

C23Sl3e** - C12C13S23 -Ci3C23e*'^ - C12S13S23 

The theoretical model H33(l is supplemented by the experimental data supplied 
by experimenters. In the quark sector one measures two kinds of parameters: 
the moduli of the unitary matrix H33|) . see |18) . under the form of a positive 
entries matrix, written with the physicists notation 




V 



ud us ub 
cd ' cs ' cb 

V v,l V? 



yid yi yi I (34) 

t1 

tb 



where u, s, b, etc. are names for quarks, and the angles of the so called standard 
unitarity triangle, j^O], denoted by a/(j)i, P/(j)2,'y/4>3, 02]. In this paper V 
denotes either a double stochastic matrix, or a set of numbers affected by errors, 
when Vij are measured in experiments. More generally the experimental data 
can be given in terms of some functions fkiVij), k = I, . . . , N , that depend on 
the V entries, or the theoretical parameters Sij and S. 

Similarly to the general case treated in the previous section, we define our 
phenomenological model as a relationship between the theoretical object H33fl 
and the experimental data (|34|l . It is given by the double stochasticity relations 
(jll|l . which now take the form 



^ V,2-l = 0, j=u,c,t 

i—d,s,b 

^.'-1 = 0' J=d,''^,b (35) 

i—u,c^t 

and by the embedding relation of a unitary matrix into the double stochastic 
set 

V = \U\^ 

that leads to the following relations 



T/2 _ „2 T/2 _ 2 2 ^r2 _ 2 2 

T/2 _ „2 2 T/2 _ 2 2 

•'cd ~ '^12''23' ^td — *12''23' 

Vcs = Ci2C?3C23 + S13S23 + 2C12C13C23S13S23 cos 6, 

= C12C23S13 + C13S23 - 2C12C13C23S13S23 cos (5, (36) 

ytl = C23S13 + Ci2C?3S23 - 2C12C13C23S13S23 cos 6, 

Vfl = C^3C23 + Ci2S?3S23 + 2C12C13C23S13S23 COS S 
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The above phenomcnological model was introduced in jl^- The relations 1)36(1 
depend only on cos (5 which has the consequence that we can restrict S to the 
interval [0, tt], this property being equivalent to the CKM matrix invariance 
under the complex conjugation, as it was shown in section 2. 

Especially for physicists we want to make a few remarks. First we stress that 
in any phenomcnological analysis one works with two distinct objects: the first 
is the theoretical one, that in our case coincides with the matrix U, Ea. (|33|l . 
which is assumed and built as a unitary matrix; the second object is provided 
by the experimental data, V = (V^), Eg . H34|l . The aim of any phenomcnologi- 
cal analysis is twofold: a) checking the consistency of data with the theoretical 
model, and, b) determination of parameters entering the theoretical model from 
the experimental data, if these ones are consistent with it. That is the rea- 
son for making a clear distinction between the theoretical quantities and the 
experimental ones, by using different symbols for denoting them. The second 
remark concerns the double stochasticity relations, Eas. ((35|l . that are consid- 
ered by (many) high energy physicists as testing the unitarity, statement which 
is wrong, since it is well known that for n > 3 there exist double stochastic 
matrices which are not unistochastic, [2]. Checking the consistency of the data 
with the theoretical model means checking the consistency of relations (|36|l . i.e. 
we have to see if the solutions of Eqs. 1)36(1 lead to physical values for the mixing 
parameters 9ij and S; and only in this case Eos. ((35(1 together with Eos. ((36(1 prove 
the unitarity property of the data. 

Let us assume for a moment that the relations ((35(1 are exactly satisfied. 
Then it is an easy matter to find from the first five relations 1(36(1 three inde- 
pendent ones which give a unique solution for the Cy , = 12, 13, 23. In other 
words, if the experimental numbers satisfy the relations 

+ + Kd = 1 

we get always a solution for Cij that is unique and depends on the three chosen 
independent parameters. Substituting this solution in the last equations one 
gets four equations for cos 5, that lead to a unique solution for it. But nobody 
guarantees us that the solution will satisfy the physical constraint 

-l<cos(5<l (37) 

The last relation gives the necessary and sufficient condition the data have to 
satisfy in order the 3x3 matrix ((34(1 comes from a unitary matrix, i.e. it is the 
consistency condition between the data and the theoretical model. 

To better understand the above considerations and see the power of the found 
criterion ((37(1 and how it works, we will give a few numerical examples, and for 
that we will use moduli entering the first two rows. We make the following 
notation: 

Vud = a, Vus = b, Vub = c, Vcd = d, Vcs = e, and Vet = f 

First we choose as independent moduli a, b, d, and e and with them form the 
square root of a double stochastic matrix 



(a b Vl-a^- b^ 



Si 



d e Vl - rf2 _ e2 (38) 
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i.e. Sf is an exact doubly stochastic matrix, where the square is taken entry 
wise, by using the Hadamard product from Unear algebra. From the relations 
(I36|l we get the solution 



ci2 = Vud = a, ci3 



C23 



(39) 



COS 5i = ^ — ' ■ : \ ' = — ^ — (40) 

In the second case we take b, c, d, and, / as independent moduli, and get 

\ 

(41) 



So 



c 

f 



V Vfo2 + C2 - d2 y/rf2 + /2 _ 52 ^1 _ c2 - /2 / 



C12 



= Vl - 52 - c2. 



Cl3 



F, C23 



VP 



(42) 



COS. 5-2 



62(62 + e2) - d2(j,2 _ ^2 ^ ^2(^2 ^ ^2)) _ ^2(^2 ^ ^2)2 

26cdVl - 62 - c2V62 + c2 - rf2 



If the data are the entries of the following double stochastic matrix 



V = 



1 2 

4 5 

_5_ J_ 

12 10 



20 
29 
60 



we get from the equations (EHJ-ljini), and 



C12 



x/3 



V3 V6 , , 4^/15 

, Ci3 = — , C23 = — , cosdi = cosd2 = ^- 
2 4 25 



(43) 



(44) 



(45) 



and the results show that the data are compatible to the existence of a unitary 
matrix. We remark that no matter how the independent moduli are chosen, Cij 
and cos (5 takes the same value, and this is a consequence of the fact that the 
properties of a double stochastic matrix do not change by multiplying it at left 
and/or right by permutation matrices. From a mathematical point of view the 
story ends here, because we can easily reconstruct the unitary matrix whose 
moduli arc given in (|44() . by using the results H45() in the unitary matrix (|33|l . 
We get 
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Hence the reconstruction algorithm of unitary matrices from the double 
stochastic ones is the following: Start with a double stochastic matrix as H44() 
and solve the system of equations (|36|l . // the numerical value for cos S satisfies 
the inequalities (|37|l . then with the values for Cij and cosS go to and find 
the corresponding unitary matrix. 

For experimental data as those recommended in ^S] the situation changes. 
For example, by using the numbers: a = 0.9738 ± 0.0005, b = 0.22 ± 0.0026, c = 
0.00367 ± 0.00047, d = 0.224 ± 0.012, e = 0.996 ± 0.013, / = 0.0423 ± 0.0015, to 
define two doubly stochastic matrices Si and 5*2 one gets 

cos(5+ = -0.03i, cos<5J = 1.59, cos (5i~ = 1.08 (47) 
cos 5+ = 8.95 i, cos 5^ = 5.985 i, cos(5^= 7.699 (48) 

where the indexes +, c, — denote the cos S values obtained from central values+lcr, 
the central values, and, respectively, central values— Icr. The above results show 
that our criterion H37|l is very sensitive to small variations of the parameters of 
the order of errors, and, on the other hand, one sees that the fulfillment of unitar- 
ity for experimental data is not an easy problem. Wc remark that cos Si ^ cos 62, 
although both the matrices H41() and H38() lead to double stochastic matrices, but 
these ones are different because numerically, e.g. e 7^ -^/l — rf^ — /2. Hence in 
the case of experimental data we have to take care and try to find how the 
necessary and suflacient conditions for the existence of a unitary matrix could 
be implemented. In this case also the relations (|35|l arc not exactly fulfilled. 
Consequently the numbers Cij obtained from the relations (|39|l and H42() could 
be different, depending on the independent parameters we use for their deter- 
mination. On the other hand the last four relations H36|) provide us formulas 
for cos (5 and these formulas have to give the same number when comparing 
theory with experiment, by supposing the data come from a unitary matrix. 
Their explicit form depends on the independent four parameters we choose to 
parameterize the data. In fact there are 58 independent groups of four inde- 
pendent moduli that lead to 165 different expressions for cos i5. Depending on 
the explicit choice of the four independent parameters we get one, two, three or 
four different expressions for cos S. 

Looking at equations 1)3911 - 1)40(1 and 1)4211 - 1)43(1 we see that the expressions 
defining the mixing angles and phase 5 are quite different. Thus if the data 
are compatible to the existence of a unitary matrix these angles and phases 
(5(*) have to be equal, and these are the most general necessary conditions for 
unitarity; they can be written as 

< c*"^ < 1, c|"^ = c[]\ m, n 1, • • • , 58, cosJ^'^ = cosS'^^\ j = 1, . . . , 165 

The above relations are also satisfied by the double stochastic matrices, and the 
condition that separates the unitary matrices from the double stochastic ones 
is given by the relation 1(37(1 . i.e. — 1 < cos (5'*^ < 1. 

4.2 Unitarity triangle method 

The second phcnomcnological model is defined by the orthogonality relations 
of rows, and, respectively, columns of a unitary matrix, and by the double 
stochastic relations ((35(1 . Although there are six such relations, see Eas. ((12(l . 
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usually one considers only the orthogonality of the first and the third columns 
of U, relation that is written as 



Uudu:, + Ucdu:, + utdUt\ - 



(49) 



The above equation can be visualized as a triangle in the complex plane. Usually 
(|49|l is scaled by dividing it through the middle term such that the length of one 
side is 1. Taking into account that our parametrization, Eq.|j22I), of a unitary 
matrix has the entries of the first column and the first row positive quantities, 
we divide by the first term, UudU*f^^ which is positive. In fact what matters are 
the angles of the triangle, and our choice has the advantage that the triangle 
generated by Eq. (|49|1 has two angles which numerically coincide with the phases 
of U23 and C/33; together with the phases of U22 and U32 they can be used for 
the determination of the unitary matrix U, since all these angles are measurable 
quantities in experiments, see e.g. [201, or 
The other sides have the lengths 



R 



R 



(1) _ 


UcdU*, 


db.c 


UudUub' 


.(1) 


UtdU;, 


'db,t ~ 


UudU*i, 



a \/l — I 



&2 



(50) 



rf2 Va^ + 52 + cf2 _|_ g2 



&2 



where we have written on the right hand side the R-values in our choice of the 
four independent parameters by using the matrix The physical condition 

takes the form 



\R 



(1) 

db.c 



R 



(1) 

'db,t 



<l< R 



(1) 

db.c 



R 



(1) 

db.t 



(51) 



that says that with the lengths 1, Rdb,c and Rdb,t one can construct a triangle. 
If we use the matrix ()41|l we find 



R 



(2) 
'db.c 



UcdU 



cb 



R 



(2) 
db.t 



UudU*^ 
UtdU:, 



df 



:Vl-62-c2 



(52) 



UudU*^ 



P V62 +C^-(P 



■Vl-62-c2 



We remark that in Eas. l|5U|l - H52|l the left side is the same, and only the right 
side differs, because the four independent moduli we use are different. Since 
there are 58 different groups of independent moduli there will be 58 different 



,(0 



c, t, and i = 1, 



,58. 



expressions i?^^' ^ , j 

If now we use the orthogonality between the first and the second columns, 

i.e. 



UudU*u, + UcdK, + UtdUl = (53) 

one gets similarly 



(1) _ 


UcdU*, 


de 




ds^c 


UudU*, 


a b 




.(1) _ 


UtdU;, 




- a2 - d2 VI _ 52 _ g2 


'ds.t 


UudU*, 




a b 



(54) 
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and respectively 



R 



(2) 



UcdU* 



R 



(2) 
ds.t 



UudU* 



(55) 



62 



6 Vl - &2 ^ C2 



If as in the preceding case we compute the expressions on the right hand side 
using the data H44|l we find 



R 



(1) 

dh.c 



R 



db., 



7:\ R 



(1) 

db,t 



(2) 
db,t 



and respectively 



R 



(1) 

ds^c 



R 



(2) 

'ds.c 



R 



(1) 

ds.t 



R 



(2) 
ds,t 



that both satisfy the inequalities of the form (|51() . 
With the central values from one gets 



R 



R 



(2) 
db.c 



2.58 R 



(2) 
db.t 



11.72 i i?' 



(1) 

(2) 

'ds,c 



1.04 i?, 
1.016 R 



(1) 

'ds^t 
(2) 



(56) 



(57) 



0.037 i 
0.012 i 



(58) 



result that sends the same signal of incompatibility as in the previous subsection, 
see Eqs.lgHJ- 

From the above equations we can obtain the angles of the triangles generated 
by the relations H49(l and, respectively, (|53|l . For each triangle we denote by (^3 
the angle of the triangle associated to the vertex (0,0), the other two, ipi and 
(^2, being associated respectively to the vertexes (p, ??) and (1,0), where (p, 77) 
are the coordinates of the triangle apex. For the second triangle we make the 
substitution ipi —f tpi, i = 1, 2, 3, and find 



cosLpi = A\j w 0.8, cos(y92 
and respectively 

cos ■01 = ~ 



61 



10V58 



0.88, COS<y93 = --^ 
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0.19, COS ■02 = 7^7: ~ 0.65, cost's = 

20 8 



— w -0.42 (59) 



- w 0.87 (60) 



Similarly to the preceding case the necessary and sufficient conditions for 
unitarity are the constraints: all i?^-'-' > 0, and 



R 



R 



db.c 



R 



R 



db.c^ 



R 



R 



■db,t 

(«) 
ds.t 



Rdb^v *,J = 1,--- ,58 

0) 

ds.t' 



R 



,58 



(61) 
(62) 



\R 



db.c 



\R 



ds.t 



R 



R 



U) I 

■db,c\ 

U) 
ds,t\ 



<l< R 



<1< R 



db.c 

(») 

ds,t 



db.c J 

■ U) 
ds.t' 



J 



,58 
,58 



(63) 
(64) 
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where in the last equations we have written only the conditions implied by 
two orthogonality relations, although for applications we must calculate the 
constraints for all the six orthogonality relations. 

The first remark is that this approach, in the variant used by physicists, 
does not make use of the double stochasticity relations, Eas. (|35|l . the physicists 
implicitly assuming that they are satisfied by the numbers obtained from exper- 
iments. With the above notation, (p, ry), for the apex of the triangle, Eas. (|^ 
are written under the form, see e.g. [^U], |21|-|30j 



R 



R 



(1) _ 


UcdU^b 


db,c 


UudUub 


(1) 


utdu:. 


db,t ^ 


UudUub 



= v/(i-p)2 + ?f 



(65) 



The second remark is that in contradistinction to what physicists believe, p 
and r] have no special relationship with a parametrization of the CKM matrix, in 
particular that provided by Wolfenstein Indeed on the right hand, Eas. (|65|l . 
are the lengths of two sides of the above defined triangle, those ^ 1. Physical 
meaning have only the angles of that triangle which can be measured [20] , . 
The third remark is that in this approach there is no relationship between the 
CP-violating phase 6 and the angles of the unitarity triangles. Because the 
phase is interesting from a physical point of view, the phcnomenologists make 
the identification 



S = LPs, or (5 « LPs (66) 

see e.g. [21]- EDI- Looking at the numerical values obtained for cosiS = ^vTS, 
relation H45|l . and for the angles of the unitarity triangles, H59() - (|60|l . computed 
by using an exact unitary matrix, we arrive at the conclusion that the claim 
(|66|l is definitely wrong. 

A simpler example is the following: take all V^j ~ 1/3. Then all the six 
triangles are equilateral and by consequence we have (/>i = ^2 = 03 = 60°, and 
from the first phcnomcnological model we get S = 90°. 

Although this model in the form l)65|l is currently used in many phenomeno- 
logical analyzes, see e.g. j25j and |26j . it cannot provide directly numbers for 
the parameters Cy and S such that there is no (reliable) recovery algorithm for 
unitary matrices from double stochastic matrices. The positive thing is that, if 
properly used, this phcnomcnological model allows the determination of all the 
angles of all the six unitarity triangles, which are measurable quantities. Hence 
the real problem is to find a recovery algorithm for unitary matrices starting 
with measured values for all these angles. This problem was first raised by 
Aleksan et al, and in the next section we will solve it. 



5 Equivalence of the two approaches 

The relations (|48() and H58() . as well as H45() and (|59|l . have shown that the 
unitarity sends the same signal of (in)consistency between the data and the 
theoretical model, although each one in a specific way. This is natural since both 
the models are based on the unitarity property of matrices modeling the CP- 
violation. In the following we will prove that from a theoretical point of view 



19 



the above phenomcnological models are only partially equivalent in the following 
sense: if we start with four independent angles we can reconstruct more than a 
unitary matrix, the multiplicity being equal to five, and the solutions given by 
the two models are the same if and only if four moduli take the same values. This 
opens the possibility to define new phenomcnological models, in terms of moduli 
and phases, by taking into account all the experimental data. In particular we 
provide an expression for the phase i5 in the second phenomcnological model. 

The starting point is the relation (|50|l that we write in a complex form. For 
that we define the four independent angles that enter the CKM matrix (j^SJ, 
namely Uij = je'^^'J , i, j = 2, 3. In our notation, ((SOJ, UJ23 is the angle located 
at (0,0) between the positive x-axis and the complex vector U*^Ucb/U*^Uub 
oriented towards the apex (p, rj); it — UJ33 is the angle located at (1, 0) between 
the positive x-axis and the complex vector U^JJtb/U*JJub oriented also to the 
same apex. With this notation the complex form of (|50|l is written as 

Rdb,c (cos W23 + i sinw23) = tt^ (6"^) 

'^ub 

Rdb.t {- COSUJ33 + isinujss) ^ ^Jf (68) 

U„A Uub 



where from we get 



C13C23S23 smd 

smu;23 = 5 (69) 

C12 Si3 n.db,c 

C23(C12 C23 Si3 - Ci3 S23 COS 5) 

COSW23 = 5 (70) 

C12 Si3 tidb.c 

The above relations are equivalent with 

ci3 S23 sin 5 

tanw23 ; 7 (71) 

-C12 C23 si3 + ci3 S23 cos 6 

The last formula depends only on theoretical parameters entering (|33|l . and 
does not depend on the lengths of the unitarity triangle. If in (|71|l we sub- 
stitute values for Cy and cos 5 taken for example from the phenomcnological 
model, Eas. H39|) - (|40|l . we get a formula for tanw23 in terms of moduli that are 
measurable quantities. 

Conversely, from the relations H69|l - H70() we get a formula for cos 5. Indeed, 
from the identity sin^ 0^23 + cos^ 0^23 = 1 we find 

J, C23 -I- Ci3 C23 S23 — C12 Si3 R^b,c /-n\ 

COS 5 = 2 ( -^2) 

2ci2 Ci3 C^3 Si3 S23 

If in it we substitute the formulas H39|l and Rdb.c from Eas. ljSOl) we find Ea. H40fl 
for cos (5. If instead of H50|) wc use the corresponding form for Rdb.c, that comes 
from the relation H52() . one gets Ea. (|43() . and so on. The above formula shows 
that if we want to use the lengths of the unitarity triangle to obtain cos 5 we 
have to provide values for from an other source. 
In the same way from the relation (|68|1 one gets 

ci3C23sin(5 

tanw33 ■ (73) 

C12S13 S23 + C13 C23 cosd 
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and 



2 2 d2 _ 2 2 2 „ 2 2 4 
'-12*13^d6,i "-IS ^-23 ''23 '^12 *13 ''23 

2 Ci2 CIS C23 Sl3 



(74) 



The third angle is given by tt — ^23 — Of course from the relations l|H7|l - l|HH| 
we find 

— — (-cos(a;23 +W33) + » sin(w23 + W33)) = . (75) 

^db,c Ucb 

and from it we can obtain a similar formula for tan(tt'23 + W33). 
Similarly the complex form of the second triangle, Eas. (|54|l . is 



U* U 

Rds,c{cOSUJ22 +isinLJ22) = Tj*\j' C^S) 

^ud s 

u* u 

Rds,t {- COS uj32+i sin UJ32) ^ -j-rT-j7^ (77) 



where from we get 



si3S23sin^ 

tana)22 = ? (78) 

C12 C13 C23 + S13 S23 cosd 



, C23Si3Sm(5 

tana'32 = ■ (79) 

-C12 ci3 S23 + C23 si3 cos d 

From the above calculations one sees that an orthogonality relation deter- 
mines only two independent angles, and from them one cannot reconstruct the 
unitary matrix because each angle depends on the four independent parameters 
entering the generic form of a 3 x 3 unitary matrix. 

From a mathematical point of view the angles ujij, i,j = 2, 3 are not very in- 
teresting, although their existence was the essential ingredient for obtaining the 
necessary and sufficient conditions for the existence of a unitary matrix from the 
entries of a double stochastic matrix, see P] ■ From a physical point of view they 
are very interesting because they are measurable quantities. Similarly to the 
preceding cases we have to use all the six orthogonality relations, although for a 
double stochastic matrix all the angles Wy, i,j = 2, 3 have the same numerical 
values, irrespective of the orthogonality property we use. However irrespective 
what triangles we use we get the same functions tan 0^22, tana;23, tanti;32, and 
tantj33. If we use the orthogonality of the second and the third columns, respec- 
tively of the second and the third rows, we get the angles of the corresponding 
triangles as linear functions of 0^22, ^23, W32, W33, see |44) . 

In the following we give the necessary and sufficient conditions for recovery 
of a unitary matrix when we know the angles , i, j = 2,3, solving the problem 
first raised by Aleksan et al In the following we make the notation 

tanu)22 = ^22, tanw23 = t23, tancj32=t32, tancJ33=t33 (80) 
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and from the Eqs.(|7Tl),lE2Il, and ijTHll-lEni) we get 

2 ^23 ^33(^22 — ^32) 



"13 



-23 



^2 
'-12 



^22 ^23(^33 ^^32)+ ^32 ^33(^22 — ^23) 

^32 ^33(^22 — ^23) 
^22 ^23(^33 — ^32) + i32 ^33(^22 — ^23) 
iVi 
N2 



(81) 



where 

-^1 = (^22 ~ i23)(^22 ~ i32)(i23 — i33)(^32 ~ ^33) 

N2 = ^23 ^32 + ^22 ^33 + *23 ^32(^22 + ^33) + ^22 ^33(^23 + ^32) ~ (82) 
2 t22 ^23 ^32 ^33 [1 + (^23 + ^32) (^22 +^33) — ^23 ^32 ~ ^22 ^33] 

Substituting the values for Cy from relations H81|l and (|82() in any equation 
()71|1 . (|73|) . H78|I . (|79|1 . or a combination of them, we get a formula for cos (5, that is 
too long to be written down here. Hence the necessary and sufficient conditions 
for the existence of a unitary matrix coming from the angles ujij are 



< < 1, < < 1, < C23 < 1, -l<cos(5<l 



(83) 



By using the numerical values 



1 /77 _ /33 _ V23T 3V23T 

^22 — g Y ^23 — y -y, tz2 — ~ 61 

obtained from the matrix H46|) we get by using the relations H81|l and (|82|l 

C12 = f^lS = — , '^23 = — (85) 

showing that Eas. H81|l - H82() uniquely define the parameters , in perfect accord 
with 

If in Eqs.(|ni),IZSll, and (jTHJ-IZni) we substitute the angles as given by 
and H82|) we get equations for cos 5 that lead to the solutions 

cos5 - ^vl' ^°'^="iS/!' ^^^-^^^vl' ^^^^ 

cos^ = -\^|\. cosJ = -|y| (87) 

In fact from each equation H71|l - H73|) and (|78|I - H79|I one gets two solutions for 
cos (5, and only one of them coincides with that found in the first phenomeno- 
logical model, see H45|l . Hence the problem of recovering a unitary matrix when 
we know four independent angles is not unique, and the finite multiplicity 
is at least five. However this result does not contradict the general result stated 
in Theorem 2, theorem which gives oneness then and only then when we use 
four independent moduli. To see what happens in the above case, we recover 
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the unitary matrix by using Cy taken from the relations H85|) . and the second 
value for cos S 

139 /3 

cos = -1 / - 

116 V 5 

The moduli matrix has the form 







1 


1 






( ^ 


2 


6 


\ 




1 

4 


47 
1856 


1345 
1856 






, A. 

\ 12 


881 
1856 


605 
5568 





By comparing the matrix Vi with the original one, Ea. H44(l . we see that the 
elements of the first row and column coincide, the others are different. In order 
to obtain a unique solution we can use the relation ifT^ , or any other equivalent 
to, which make use of one of the lengths of unitarity triangle sides. Hence by 
using information coming only from triangles angles we have a finite multiplicity 
solution. The unicity is obtained then and only then when the information is 
supplemented by an independent modulus, e.g. in the above case V{2,2), or a 
length of a unitarity triangle. The phenomenological implications of the above 
results on the global fit methods for recovering a unitary matrix from moduli 
and angles will be treated elsewhere. 



6 Recovery of unitary matrices from experimen- 
tal data 

If the data come from an exact numerical matrix the problem to solve is quite 
simple: we have to test the stochasticity property, and, afterwords, the unitarity 
constraints, i.e. the condition, —1 < cos (5 < 1, in the unitarity condition 
method, or the inequalities, |i?a/3,7 — Raf3,-(\ < 1 < Rap.-y + Rai3,-y, coming 
from two orthogonality relations in the case of the standard unitarity triangles 
approach. If the data pass anyone of the tests, one can easily reconstruct the 
unitary matrix Ea. H33|l from the data (|34|) . as the numerical examples from 
the previous sections have shown. If the physical conditions are violated, there 
is no compatibility and the discussion ends here. The real problem, from a 
physical point of view, is when the data come from experiment, i.e. are numbers 
affected by errors, and the problem is how we proceed in this situation, because 
neither the double double stochasticity relations, nor the unitarity constraints 
are exactly satisfied. 

There is the place where the gauge invariance subgroup K of unitary matrices 
enters the game, and its implications are the following. We have to find all the 
four independent moduli groups and find all the possible forms for cos 5, and, 
respectively, for the lengths of the unitarity triangles. And we have to impose 
that the numerical values for them should be approximately equal. The usual 
case with the present data is that the numerical values obtained for moduli 
are such that cosf^^'^ ^ cosd^^\ and/or i?^^ ^ ^ R-ap -y^ * J- Even more 
cos 5^*^ could be outside the physical region, or the lengths of unitarity triangles 
are imaginary, or if they are real are not compatible with the existence of a 
triangle, as the numerical examples provided in the paper show. Hence, in 
contradistinction with the nowadays usage, see |2U) . we have to devise a fitting 
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model that should implement the fulfillment of the above theoretical constraints, 
and which should take into account the experimental data. 

The method we expose here is discussed in more detail in It is a least 
squares method for checking the compatibility of the data with the theoretical 
models in both the approaches, and if the data pass the physical conditions 
imposed by unitarity, from the fits one gets values for the parameters entering 
the theoretical model that by assumption gives a reliable description of the 
physical reality that is investigated by experiment. 

It follows that in both the approaches the x^-function must contain two sep- 
arate terms, the first have to impose the fulfillment of the unitarity constraints 
by the free parameters entering the physical model, and their best determina- 
tion from the data, and the second one should depend on physical quantities 
that are measured in (different) experiments. Thus our proposal for the first 
terms is 



i<j j—u,c,t \i—d.,s,b j 

+ E ( E ^^'-ij ' -l<cos<S«<l (89) 

j—d,s,b \i—u,c,t I 

for the unitarity condition method, and, respectively, 



= E(<.-<.)^+ E E + E E 

j—u,c,t \ i—d,s,b J j—d,s,b \i—u,c,t 

i<j 

-"q/3,7 - l-"'Q/3,7 ^a/3,-(\ -"a/3,7 ^ a/3,7 

for the unitarity triangles method. Both the Xi 2 formulas test the double 
stochasticity property and the unitarity; from the point of view of numerical 
computation the unitarity property is the most difficult to satisfy. 
Concerning the second component of x^-tcst it is of the form 



(91) 



where di are theoretical functions depending on the theoretical parameters Sij 
and S, or on Vki, or/and the angles which are given by the phenomeno- 

logical model one works with, while di are the measured experimental data for 
di, and a is the vector of errors associated to di. The formulas 

2 2,2 

and, respectively, 

2 2,2 
Xt =X2+ X3 
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are our proposals for x^-tools necessary in analyzing the experimental data. 

A remark is the following: as we have seen the second phcnomenological 
model does not provide a formula for cos 6. However any global fit done by using 
either Xu Xt gives values for all the moduli. Hence for the reconstruction of a 
unitary matrix in the second phcnomenological model we have to use formulas 
such as (I74II . A true global fit will be that which will use all the experimental 
data by merging the above two phenomenological models. 

The convexity property of the double stochastic matrices allows us to devise 
a method for doing statistics on unitary matrices that is still an open problem in 
the physical literature. Let suppose that by doing a fit with the above methods 
we got n moduli matrices that are consistent with n (approximate) unitary 
matrices, Ui, U2 - ■ ■ ,Un- The convexity property together with the embedding 
(0) tell us that the matrix 

n n 

1=1 1=1 

is double stochastic, where we use the Hadamard product, so the above relation 
is understood as working entry wise. Then the correct formulas for the mean 
value < M >, and the error matrix ctm are given by 



i=l 



(E m')/n- < M >4 (92) 

If the entries of the mean value matrix, < A/ >, obtained in this way are not 
too far from the entries coming from a unitary matrix, one can reconstruct from 
< M > an (approximate) unitary matrix by using the technique developed in 
the paper. 

7 Conclusions 

Our main reason for studying the separation criteria between the double stochas- 
tic matrices and the unistochastic ones was their importance in high energy 
physics, where the nowadays algorithms for reconstruction of unitary matrices 
from experimental data are not yet reliable in our opinion. Fortunately for the 
3x3 matrices, that seems to be the physical choice in the electroweak inter- 
action, there are explicit formulas for the independent parameters entering a 
unitary matrix in terms of four independent elements of a double stochastic 
one. That allows us to check the unitarity properties of exact double stochastic 
matrices, and an easy reconstruction of the unitary one from the entries of the 
double stochastic matrix when the compatibility conditions are fulfilled. These 
formulas can be used to define x^-functions for checking the compatibility be- 
tween the experimental data and the unitarity property of the CKM matrix, and 
to recover a unitary matrix from error affected data. More important, starting 
from the convexity of the Birkhoff's polytopc, we found a method for doing 
statistics on the (moduli of) unitary matrices. 



<M> 



\ 



\ 
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We have also shown that, because the unitarity triangles method [201 and 
the unitarity condition method |43| . being both consequences of the unitarity 
property, are completely equivalent when and only when they are formulated in 
terms of four independent moduli. In the same time wc have shown that the 
unitarity triangles method has to make effective use of the double stochasticity 
relations in order to obtain reliable results. Writing the unitarity triangles 
method in complex form wc have obtained formulas for the four independent 
phases entering a unitary matrix, these phases being the angles of the unitarity 
triangles. This opens the possibility to treat coherently all the experimental 
data available on moduli and angles by merging the above phenomenological 
models into a true global one, the aim in view being a precise determination of 
the phase S that is the key parameter in understanding the CP-violation. 
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